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Managementuittreksel 

Titei  Speckle  reduction  in  SAR  Imagery  by 

various  multi-look  techniques 
Auteur(s)  Ir.  F.P.Ph.  de  Vries 

Datum  :  januari  1998 

Opdrachtnr.  :  A93KM727 

IWP-nr.  767.2 

Rapportnr.  :  FEL-96-A015 


SAR  beelden  bevatten  speckle,  een  type  ruis  dat  ontstaat  door  de  coherente 
verwerking  van  radar-echo’s  van  gedistribueerde  doelen.  De  aanwezigheid  van 
deze  ruis  vermindert  de  detecteerbaarheid  van  objecten  in  het  beeld.  De  mogelijk- 
heid  om  gedistribueerde  doelen  te  scheiden  en  te  classificeren,  vermindert  even- 
eens.  Zoals  elke  ruisvorm  kan  speckle  worden  gereduceerd  door  middelen,  het 
samennemen  van  een  aantal  beeldrealisaties  van  dezelfde  scene.  Een  van  de  tech- 
nieken  die  kan  worden  gebruikt  is  ‘multi-look  processing’.  Hierbij  worden  de 
samenstellende  beelden,  de  ‘single-look’  beelden,  gevormd  door  voor  elk  beeld 
een  deel  van  de  beschikbare  signaalbandbreedte  te  gebruiken.  Elk  deel  wordt  dan 
verwerkt  tot  een  single-look  beeld  van  dezelfde  scene.  Het  ruisgereduceerde  beeld, 
het  multi-look  beeld,  wordt  verkregen  door  de  beelden  incoherent,  pixel  voor 
pixel,  op  te  tellen. 

Het  speckle-gereduceerde  beeld  heeft  nu  een  verbeterde  radiometrische  resolutie. 
Echter,  omdat  de  single-look  beelden  slechts  een  deel  van  de  beschikbare  band- 
breedte  gebruiken,  zal  de  geometrische  resolutie  slechter  zijn.  De  verbeterde  radio¬ 
metrische  resolutie  wordt  verkregen  ten  koste  van  de  geometrische  resolutie.  De 
resulterende  resoluties  worden  bepaald  door  de  manier  waarop  de  beschikbare 
bandbreedte  wordt  verdeeld  en  de  gebruikte  weegfunctie.  Het  toepassen  van  weeg- 
functies  verbetert  het  zijlusgedrag  van  de  corresponderende  impulsresponsie,  de 
lekkage,  het  zich  uiten  van  (sterke)  reflecties  buiten  de  doelpositie  wordt  dan 
minder.  Dit  verbeterde  ‘ lekkage’ gedrag  gaat  ten  koste  van  de  geometrische 
resolutie,  dat  wil  zeggen  de  breedte  van  de  hoofdlus  van  de  puntdoelresponsie 
neemt  toe. 

In  dit  rapport  worden  een  aantal  multi-look  technieken  onderzocht,  die  zich  onder- 
scheiden  in  de  vorm  van  de  weegfunctie,  de  breedte  van  de  weegfunctie,  de 
combinatie  van  verschillende  weegfuncties  en  de  posities  van  de  weegfuncties  ten 
opzichte  van  elkaar. 

De  voor  het  beoordelen  van  de  technieken  benodigde  maten  worden  behandeld. 
Deze  maten  zijn:  het  equivalente  aantal  beelden,  ENL,  hetgeen  een  maat  is  voor  de 
ruisvermindering,  en  de  vermogenshalfwaardebreedte  van  de  hoofdlus  van  de 
puntdoelresponsie,  zijnde  een  maat  voor  de  geometrische  resolutie. 
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1.  Introduction 


Speckle  appears  in  S  AR  images  due  to  coherent  processing  of  radar  echoes  of 
distributed  target  scatterers.  Its  presence  in  SAR  images  reduces  the  detectability 
of  objects  in  the  image  and  also  the  capability  to  separate  and  classify  distributed 
targets.  Noise  can  be  reduced  by  averaging.  One  technique  which  can  be  used  is 
multi-look  processing.  Here  the  available  signal  bandwidth  is  divided  in  L  parts. 
Each  part  is  processed  yielding  a  single  look  image  of  the  same  scene. 

By  adding  the  looks  incoherently  pixel  by  pixel  the  speckle  reduced  image  is 
obtained.  The  speckle  reduced  image  has  an  improved  radiometric  resolution. 
However,  since  the  single  look  images  use  only  part  of  the  available  bandwidth, 
the  geometric  resolution  will  be  worse.  Improvement  of  the  radiometric  resolution 
is  obtained  at  the  cost  of  the  geometric  resolution.  The  resulting  resolutions  are 
determined  by:  how  the  available  bandwidth  is  partioned  and  the  weighting 
function  used.  Using  a  weighting  function  improves  the  sidelobe  behaviour  of  the 
corresponding  impulse  response  at  the  expense  of  a  broadening  of  the  main  lobe 
and  thus  a  degradation  of  the  geometric  resolution,  see  [3]. 

In  section  2  we  introduce  the  improvement  factor  ENL,  the  Equivalent  Number  of 
Looks.  The  various  multi-looking  techniques  are  treated  in  section  3. 
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2.  The  improvement  factor  ENL,  the  Equivalent 
Number  of  Looks 


The  SAR  image  of  an  area  of  distributed  targets  will  show  speckle  noise.  The 
observed  intensity  from  a  pixel  follows  exponential  statistics.  The  exponential 
distribution  of  a  single  look  intensity  image  is  given  by 

-I 

p(I)=-^e2<^^ 

2a 

where  the  mean  value 
mi  =  E[I]  =  2a2 


and  the  standard  deviation  is  equal  to  this  mean  value 


a,  =  E[(l-I)' 


=  2a2 


In  the  sequel  we  assume  an  available  processing  bandwidth  of  Ba  with  correspond¬ 
ing  available  azimuth  integration  time  Ta* 

We  further  assume  a  flat  azimuth  SAR  spectrum,  see  Figure  1.  In  practice  this 
assumption  is  not  valid  since  the  azimuth  spectrum  is  shaped  by  the  azimuth 
antenna  gain.  In  satellite  SAR  the  available  bandwidth  Ba  is  roughly  0.8  PRF  and 
the  edges  of  the  spectral  interval  coincide  with  the  3-dB  points  of  the  azimuth 
beam.  We  also  assume  that  the  thermal  noise  is  small  compared  to  the  clutter. 


Pc(f) 


Figure  I:  The  azimuth  SAR  spectrum. 

The  noise  behaviour  of  single  look  images  is  the  same,  irrespective  of  their 
geometric  resolution  and  their  spectral  position  in  the  available  bandwidth  (flat 
spectrum  assumption!).  To  show  this,  let  the  mean  intensity  value  of  look  1  be  mi 
and  the  variance  vari  and  the  corresponding  values  of  look  2  m2  and  var2,  then  by 
a  simple  scaling  we  can  obtain  a  scaled  image  with  the  mean  value  m\. 
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m2'  =  km2  =mi 

where  k  is  the  scaling  constant.  For  the  variance  of  the  scaled  image  we  find 

var2'  =  k^  var2 

with  var2  =m2  we  have 

var2'  =  k^m2^  =  m|^  =  vari 

Thus  the  scaled  image  has  the  same  mean  value  and  variance  as  the  look  1  image. 


We  now  compare  the  variances  of  a  multi-look  image  with  m2  and  var2  composed 
of  single  look  images  with  mj  and  vari.  The  multi-look  image  is  formed  so  that 
m2  =  mi.  We  define  the  equivalent  number  of  looks  (ENL)  by 


var2 


1 

ENL 


vari 


Thus  ENL  = 


vari 


m, 


var2  var2  var2 

The  ENL  corresponds  to  the  effective  number  of  statistically  independent  looks 
and  is  given  by 


ENL: 


E[In 


var[lM] 

where  is  the  intensity  of  the  multi-look  image. 


In  the  next  chapter  we  discuss  several  multi-look  techniques  where  a  multi-look 
intensity  image  is  a  weighted  sum  of  single  look  intensity  images. 
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3,  Multi-look  techniques 

3.1  Rectangular  weighting  function,  equal  looks 

The  rectangular  weighting  function  is  given  by 
w(t)  =  l  -T/2<t<T/2 

Figure  2  shows  the  weighting  function  and  the  corresponding  azimuth  impulse 
response  profile  g(ta). 


0.8 


0.2“  . 

O' - ’ - ^ 1 - 1 _ I _ I _ I _ I _ 

-0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 


3.1.1  Full  resolution  image 

The  full  resolution  image  uses  the  whole  available  bandwidth  and  no  weighting 
function.  So  the  best  geometric  resolution  is  obtained.  Using  the  3-dB  resolution 
measure,  the  3-dB  bandwidth  of  the  azimuth  impulse  response  is  0.8  l/B^,  where 
Ba  is  the  available  bandwidth. 

Since  no  geometrical  resolution  is  sacrificed,  no  speckle  reduction  is  obtained. 
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ENL  = 


E[I] 

var[l] 


=  1 


3.1.2  Combination  of  Li  independent  looks 

The  available  integration  time  is  divided  into  Lj  discrete  contiguous  segments. 
Each  segment  is  used  to  generate  a  lower  resolution  single  look  image.  The  width 

T 

of  each  segment  is  thus  T  =  — .  The  geometrical  resolution  is  reduced  by  a  factor 

Li 

Lj,  Thus  the  3-dB  bandwidth  of  the  azimuth  impulse  response  is 

Ba 

The  mean  values  and  (cross-)variances  of  the  look  images  are  given  by 
Ii  =  I  for  each  look  image  1. 

(lj -ij)  =varl  for]  =  i  and 

=  0  for  j  i  (independent  looks!) 


The  multi-look  intensity  is  given  by 

Li 

Im  = 

1=1 

Li 

The  mean  value  Im  =  cx^^Ii  =aLiI, 

1=1 

The  mean  value  is  chosen  equal  to  the  single  look’s  mean  value,  so 
aLi  =1 


The  variance  follows  from 


2  r 

(Im  -Im)  =a2  X(l'-i') 
Vi=i 

Hence 


=  a^Li  var(I) 


variM  = 


var(I) 

Li 


Thus  the  equivalent  number  of  looks  is 


ENL  =  -^  =  Li 
variM 


Hence  the  equivalent  number  of  looks  is  equal  to  the  number  of  independent  looks. 

3.1.3  Combination  of  L  overlapping  looks 


The  available  integration  time  is  divided  into  L  overlapping  segments.  The  width 
of  each  segment  corresponds  to  one  of  the  Li  looks  in  the  independent  case  (see 
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3.1.2),  T  =  Ta/Lj.  The  L  single-look  images  are  incoherently  averaged.  The  mean 
values  and  the  (cross-)variances  of  the  look  images  are  given  by 

ii  =  I  for  each  look  image  1. 


var(Ii )  =  (li  -  Ii  =  Oi^  =  var  I  for  each  look  image  1 


(li  -Ij)  =  pyOiOj  =  pij  varl 

where  py  is  the  correlation  coefficient  between  the  i-th  and  j-th  look  image.  The 
correlation  coefficient  py  depends  on  the  degree  of  overlap. 

Since  all  looks  have  equal  weighting,  0  <  py  <l(py  =1,  complete  overlap  and 
Py  =  0,  complete  separation).  In  appendix  A  the  formula  for  the  correlation 
coefficient  is  derived.  Figure  3  shows  the  correlation  coefficient  as  a  function  of 
the  fractional  overlap.  As  is  shown  in  the  appendix,  the  correlation  coefficient 

Pij  =  piiij  between  the  intensity  images  is  equal  to  |pxiXj|^ ,  where  is  the 
correlation  coefficient  between  the  complex  valued  images  Xj  and  xj. 

The  multi-look  intensity  is  given  by 

L 

Im  =CX^Ii 
1=1 

Proceeding  as  before  we  arrive  at 

?  T  1 

Im  =  I  oc  =  — 

L 


var(lM) 


i=lH 


u 


varl 


Thus 


ENL  = 


U 

L  L 


It  turns  out  that  an  increase  of  the  ENL  is  possible,  see  figure  4. 

Since  the  width  of  the  looks  is  the  same  as  in  case  3.1.2,  the  geometric  resolution 
is  the  same. 
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Rectangular  weighting  function 


Figure  3:  The  correlation  coefficient  as  a  function  of  the  fractional  overlap 

complex  images  \  p^y  (x)  |  - 

intensity  images  pi  (x)  . . 


Equivalent  Number  of  looks,  Li  =  2  and  4,  rectangular  weighting 


3.2  Hamming  weighting  function,  equal  looks 

The  Hamming  window  is  a  frequently  used  one  (see  3).  The  function  is  given  by 
w(t)  =  a  +  (l-a)cos^-Ytj  -T/2<t<T/2 

with  a  =  0.54. 

The  3-dB  bandwidth  of  the  azimuth  impulse  response  profile  is  1 .30/Ba. 


20log10(lg(ta)l)  w(t/T) 
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Figure  5  shows  the  weighting  function  and  the  corresponding  azimuth  impulse 
response  profile  g(ta). 


t/T 


3.2.1  Combination  of  Lj  independent  looks 


Proceeding  as  in  3.1.2  we  find 


Li 


The  multi-look  intensity  image  is  Im  =  —  ,  the  mean  value,  Im  =  I  and  the 

1=1 

variance,  var(lM )  =  ,  thus  the  ENL  =  L\. 

Li 


Comparing  this  case,  the  Hamming  window,  with  the  results  from  3.1.2,  the 
rectangular  window,  we  see  that  the  speckle  reduction  is  the  same.  However,  the 
geometric  resolution  is  now  worse.  So,  we  expect  that  with  Hamming  windows  a 
better  speckle  reduction  is  attainable.  The  technique  discussed  next  shows  that  this 
is  indeed  possible. 


rhol(x)  and  Irhoxyl 
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3.2.2  Combination  of  L  overlapping  looks 

Proceeding  as  in  3.1.3  we  find 
The  multi-look  intensity  image  is 

^  1=1 

the  mean  value  Im  =I 

L  L 

^]£pijvarl 
the  variance  var(lM)=  - 

Jj 

the  Equivalent  Number  of  Looks  is  ENL  =  — — ^ - 

XXP'j 

i=i  j=i 

Figure  6  shows  the  correlation  coefficient  as  a  function  of  the  fractional  overlap. 


Hamming  weighting  function,  a  =  0.54 


intensity  images  Pi  (x) 


As  derived  in  appendix  A,  py  =  pi^i,  =|pXiXj|^,  where  pi^i.  is  the  correlation 
coefficient  of  the  intensity  images  and  Px^xj  is  the  correlation  coefficient  of  the 
complex  images. 

Comparing  the  p-values  from  this  figure  with  those  of  figure  3,  we  see  that  the 
‘Hamming’  values  are  smaller.  This  means  that  with  the  same  value  of  the  sum  of 
the  correlation  values  we  can  have  more  looks  in  the  Hamming  case,  thus  an 
increased  ENL. 
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Figure  7  shows  the  ENL  as  a  function  of  the  (overlapped)  looks.  In  the  sequel  we 
refer  to  this  multi-look  technique  as  the  traditional  technique. 


Equivalent  Number  of  looks,  Li  =  2  and  4,  Hamming  weighting,  a  =  0.54 


fractional  overlap  x 

Figure  7:  The  ENL  as  a  function  of  the  fractional  overlap. 


3.3  Hamming  weighting  function^  two  types  of  looks 

The  sections  3.1.3  and  3.2.2  are  concerned  with  overlapping  looks,  3.1.3  with 
looks  computed  with  rectangular  weighting  functions  and  3.2.2  computed  with 
Hamming  weighting  functions.  In  both  cases  an  ENL  value  is  obtained  larger  than 
the  number  of  independent  looks.  As  explained  this  originates  from  the  small 
values  of  the  correlation  coefficient  values  when  the  overlap  is  not  too  large.  Since 
in  the  Hamming  case  an  inefficient  use  of  the  signal  spectrum  is  made  by  the 
independent  looks,  a  large  number  of  looks  can  be  inserted  with  comparably  small 
p-values.  The  ‘holes’  between  the  independent  looks  can  profitably  be  filled  with 
low  p-value  looks.  However,  the  bandwidth  at  the  extremities  of  the  available 
bandwidth  are  still  not  well  covered.  In  reference  [1]  this  is  remedied  by  a 
technique  called  the  Improved  Multi-Look  technique  (IML-technique).  The 
extreme  bandwidth  is  covered  by  smaller  looks.  The  time  interval  length  of  a  small 
look  is  the  half  of  the  length  of  a  large  one.  Figure  8  shows  how  the  looks  are 
positioned  in  the  spectral  interval.  In  the  next  section,  we  derive  the  increased  ENL 
and  the  decreased  geometrical  resolution. 
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Figure  8:  Look  configuration  of  the  IML  technique,  Li  -  2,  overlap  =  50%. 

3.3.1  Ll  overlapping  large  looks  and  Ls  overlapping  small  looks 

The  mean  values  and  the  (cross-)variances  of  the  look  images  are  given  by 
iui  =iL  fori  =  1,2 . Ll 

i,.  =Is=Il/2  for  j  =  1,2, . Ls 

The  (crosS")variances  are 

(iti -iLi)(lLj =  =PiLiiLj  var(lL) 
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The  multi-look  intensity  image  is 

Ll  Ls 

Im  +0Cs^Isj 

i=i  j=i 

Thus,  the  mean  value 


Ls^ 

Im  +asX  Isj=aLXiLi  +as2^Isj  =aLLLlL+asLs^ 

i=l  j=l  i=l  j=l  ^ 

This  mean  value  is  chosen  equal  to  the  mean  value  of  a  large  look. 
Thus 

WlLlIl  +  (XsLs  —  =  II 


Hence 


cclLl  + 


c^sLs 

2 


=  1 


The  variance  is  obtained  by 


var(lM)-(lM  -  Im)  “lLi)  +  as^(lsj 


i=l 


Ll  Ll  Ll  Ls  Ls  Ls 

+2aLas^^ai,,ai5.pi^,i3.  +as^XX^‘si^isiPis,isj 

i=l  j=l  i=l  j=I  i  =  l  j=l 


=  var(lL) 


Ll  Ll  Ll  Ls  2  Ls  Ls 

^L^XP^LilLj 

i=l  j=l  i=l  j=l  ^  i=l  j=I 


Hence 


ENL  = 


LlLl  iiJlS,  fY  ^ 

«l^XXpim.4  +«,as;££p,„.,  +-f  XXp.,., 


ik  is. 


2  Ls  Ls 


i=l  j=l 


i=l  j=I 


i=l  j=l 


The  correlation  factor  p  =  |  Pxy  I  between  two  looks  is  a  function  of  the  weighting 
functions,  the  look  sizes  and  the  overlap.  If  we  use  an  overlapping  of  50%  for  both 
the  large  and  the  small  looks,  seven  values  of  p  can  occur  (see  figure  8). 

1 .  the  looks  do  not  overlap,  |  pxy  |  =0 

2.  the  looks  are  equal  and  completely  overlap  p  =  1 

3.  the  large  looks  are  correlated,  pL 

4.  the  small  looks  are  correlated,  ps 

5.  the  large  look  and  the  small  one  are  weakly  correlated  psLi  (see  figure  8, 
second  small  look  from  the  left  and  the  second  large  look  from  the  left) 

6.  the  large  look  and  the  small  one  are  correlated  psL2  (first  small  and  large  look) 

7.  the  large  look  and  the  small  one  are  strongly  correlated  psL3  (second  small  look 
and  the  first  large  look  from  the  left). 


As  explained  in  appendix  A,  py  =  pij.  =|pxjxJ^ 


ENL  3'dB  BW 
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Substituting  these  values  we  obtain 

ENL  = - i - 5 - 

®l[Ll  +2(Ll  -l)pL]  +  aLas[2(LL  -l)psu  +2LlPsl2  +LlPsl3]  +  “^[Ls  +2(Ls  -l)Ps^] 

Using  the  equations  from  appendix  A  to  calculate  the  correlation  coefficient,  we 
obtain 

Ps^  =  0.0547 ,  pI  =  0.0547 ,  Psl,  =  0.0062 ,  pslj  =  0.2693  and 
Psl3  =  0.7746 

The  number  of  overlapping  looks,  given  the  number  of  independent  looks  Lj,  is 
given  by 


LL  =  2Li-l 
Ls  =  4Li-l 


As  an  example  we  consider  the  case  for  L,  =  4. 


Improved  Multi-Look  Technique 


alphas 

Geometric  and  radiometric  resolution 


alphas 


Figure  9:  The  Equivalent  Number  of  Looks  and  the  3-dB  resolution  as  a  function  of  ag. 


TNO  report 


FEL-96-A015 


17 


As  shown  by  figure  9,  increasing  the  value  of  as  leads  indeed  to  an  increase  of  the 
ENL  and  a  decrease  of  the  3»dB  BW. 

However,  this  result  is  predictable  as  higher  as  values  correspond  with  an 
increasing  relative  contribution  of  the  small  looks,  that  is  looks  with  a  lower 
geometric  resolution  and  we  know  that  geometric  resolution  can  be  traded  for 
radiometric  resolution.  For  as  =  0  we  have  the  result  of  section  3.2.2  with  Lj  =  4 
(3-dB  BW  =1.30/Ba).  For  the  other  extreme  as  value,  corresponding  to  aL  =  0,  we 
have  the  result  of  section  3.2.2  with  Li  =  8  (3-dB  BW  =  2xl.30/Ba). 

As  a  matter  of  fact,  the  technique  offers  to  trade  radiometric  and  geometric 
resolution  in  a  continuous  way  via  the  parameter  as. 

The  technique  fills  the  ‘spectral  dips’  both  in  the  middle  part  of  the  spectral  range 
and  at  the  extremities.  The  filling  at  the  extremes  sets  this  technique  apart  from 
other  techniques.  In  order  to  prove  the  favourable  effect  of  filling  the  edges,  we 
compare  the  technique  with  one  that  fills  the  middle  part  in  an  optimum  way.  The 
latter  technique  uses  looks  each  having  the  same  size  T  (and  Hamming  weighting 
function).  The  available  bandwidth  is  completely  filled  with  L  overlapping  looks. 
The  number  of  looks  is  taken  such  that  the  maximum  ENL  occurs.  Both  techniques 
yield  the  same  geometrical  resolution;  this  is  done  by  choosing  the  corresponding 
values  of  as  and  T.  The  comparison  procedure  is 

•  chose  as  this  yields  the  ENLiml  (as)  and  3-dB  BW  (as) 

•  take  the  corresponding  3-dB  BW  (as)  value  to  determine  the  Hamming  window 
size  T(as) 

•  compute  the  optimum  L  value  for  this  window,  that  is  the  value  of  L  that  yields 
the  maximum  ENLtrad- 

The  results  are  shown  in  figure  10.  From  this  figure  we  see  that  for  almost  all 
values  as  the  IML  technique  yields  a  higher  ENL  value.  The  exceptions  are  near 
the  edges  because  there  the  techniques  are  equal  (they  use  the  same  window  size 
T),  except  for  the  overlap  which  is  optimum  for  the  traditional  technique  and  50% 
for  the  IML  technique. 

The  quotient  reaches  a  maximum  for  as  =  0.08. 

ENL  trad 

Figure  14  also  shows  that  for  the  same  geometrical  resolution  the  IML  technique 
performs  better  than  the  traditional. 

Figures  11a  and  11b  show  the  azimuth  impulse  response  of  the  IML  technique  for 
three  values  of  as.  The  value  as  =  0  yields  the  result  of  a  single  look  image  with 
Lj  =  4.  The  other  extreme  value  (aL  =  0)  yields  the  results  of  a  single  look  image 
with  Li  =  8.  The  bottom  curve  shows  the  result  for  as  =  0.08.  This  result  is 
especially  interesting  in  that  it  shows  the  relative  meaning  of  the  geometric 
resolution  measure  in  this  case.  Looking  at  the  main  lobe  we  see  that  only  the 
upper  part  has  the  character  of  a  lobe  of  a  single  look  with  the  same  3-dB  value. 
The  bottom  part  corresponds  with  the  lower  geometrical  resolution  of  the  smaller 
look. 
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Figure  JO:  Comparison  of  the  IML  and  traditional  technique  with  optimum  overlap. 


-3 

time  arbitr 


0 

3  -50 


-100^ 

-4 


-3 

time  arbitr 


-3 

time  arbitrar 


Figure  I  la:  Azimuth  impulse  re. 
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Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  15  looks  IS 


Figure  lib:  Azimuth  impulse  response  of  the  IML  technique  for  several  values  of  ag,  linear  scale. 

3.3.2  Ll  overlapping  large  looks  and  2  small  looks 

Figure  12  shows  the  configuration  of  the  technique.  Comparing  with  figure  8  from 
the  IML  technique  in  section  3.3.1,  we  see  that  the  small  looks  in  the  middle  part 
of  the  available  bandwidth  are  absent.  The  idea  is  that  the  radiometric  resolution  is 
improved  by  the  two  edge  looks  but  that  the  geometrical  resolution  deteriorates 
less  because  of  the  small  number  of  small  looks. 

The  mean  values  and  the  (cross-)variances  of  the  look  images  are  given  by 
Iti  =Il  for  i  =  1,2, . Ll 

Is,=Is4 


for  j  =  1,2 
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The  cross  variances  are 

(II,  -lLi)(lLj  =  =PiLiiLj  var(lL) 

(ILi  -lLi)(lSj  -iSj)=PlLjIs.aiLiCls.  =Pli.iIsj<JlL  “~  =  PlLilsj  - - 

(ls.-is.)(ls,-Is,)  =  Ofori;^j 

^2 

=  OisiCTis.  =^  =  var(lL)fori  =  j 

Proceeding  as  in  section  3.3.1  we  arrive  at 
Im  =  II 


o^lLl  + 


o^sLs 

2 


=  1 


and 


ENL  = 


1 


Ll  Ll 


Ll  Ll 


“l££p'^i'4  +«Las££piu.si 


i=l  j=:I 


i=l  j=l 


i-l  j=I 


As  an  example  we  consider  the  case  for  overlap  =  50%,  number  of  independent 
looks  Lj  =  4,  Ll  =  2Li-l,  Ls  =  2. 

Figure  13  shows  the  ENL  and  3-dB  BW  as  a  function  of  as,  compare  this  figure 
with  figure  9.  Figure  14  shows  the  ENL  as  a  function  of  the  3-dB  BW;  also  shown 
are  the  curves  of  the  traditional  and  IML  technique. 


Figure  12: 


Look  configuration  of  the  technique,  Li  =  4,  50%  overlap. 
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LI  =  7,  50%  overlap,  large  looks,  Ls  =  2,  small  edge  looks 


alphas 


Figure  13:  The  Equivalent  Number  of  Looks  and  the  3-dB  resolution  as  a  function  of  cxg. 

From  this  figure  it  can  be  concluded  that  this  multi-look  technique  is  not 
appreciably  better  than  the  IML  technique.  Only  for  the  lower  values  it  is  slightly 
better  and  for  the  most  part  it  is  unusable.  Note  that  in  the  3-dB  BW  range,  where 
it  performs  almost  equally  well  as  the  IML  technique,  it  does  so  in  a  far  more 
efficient  way,  only  9  looks  have  to  be  computed  compared  with  22  in  the  IML 
case. 

Figures  15a  and  15b  show  the  azimuth  impulse  response  of  the  technique  for  three 
values  of  as-  The  value  as  =  0  yields  the  result  of  a  single  look  image  with 
T  =  Ta/4.  The  other  extreme  value  (ttL  =  0)  yields  the  result  of  a  single  look  image 
with  T  =  Ta/8.  The  bottom  curve  shows  the  result  for  as  =  0.2  which  corresponds 
with  the  maximum  ENL  value  attainable.  This  curve  shows,  as  its  corresponding 
IML  curve  in  figures  1  la/b,  the  relative  meaning  of  the  geometric  resolution 
measure  in  this  case.  Observing  the  main  lobe,  we  see  that  only  the  upper  part  has 
the  character  of  a  lobe  of  a  single  look  with  the  same  3-dB  BW  value.  The  bottom 
part  corresponds  with  the  lower  geometrical  resolution  of  the  smaller  look. 
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Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 


time  arbitrary  unit  ->  aL  =  0  aS  =  1  3-dB  BW  =  2.6875  /  TL 


Figure  1 5a:  Azimuth  impulse  response  of  the  L7S2  multi-look  technique  for  several  values  of  ag,  logarithmic 
scale. 
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Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 


Figure  15b:  Azimuth  impulse  response  of  the  L7S2  multi-look  technique  for  several  values  of  ag,  linear  scale. 


3.4  Two  types  of  looks,  different  in  weighting  function  and  look 
size 

3.4.1  Ll  overlapping  large  looks  and  Ls  overlapping  small  looks 

We  proceed  along  the  same  lines  as  in  section  3.3.1.  The  only  difference  being  the 
mean  value  intensity  of  the  small  look.  In  3.3.1  we  considered  the  case  where  both 
weighting  functions  are  Hamming  ones  and  the  look  size  of  the  small  look  is  one 
half  of  that  of  the  large  look,  thus  Is  =  II  /  2 .  In  the  general  case  the  number  2  has 
to  be  replaced  by  another  number,  here  called  pm*  the  value  we  derive  below. 

The  mean  intensity  value  is 
r/2 

-Tfl 
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where 

Pc  is  the  clutter  power  spectral  density 
w(t)  is  the  weighting  function  and 
T  is  the  length  of  the  weighting  function. 

Evaluating  the  mean  values  give 
Ir  =P,Tr  andlH  =  0.375 

for  the  rectangular  and  Hamming  windows  respectively. 
Thus 


P  =  it  =  il.  =  0.3974^ 

Is  Ir  Tr 

For  the  resolution  holds 

Br  =-?^and  Bh  =  — 

Tr  “  Th 

k  -  OQQ 

with  —  =  -^ — ,  see  ref  [3]  table  I. 
kH  1.30 

k 

Requiring  equal  resolution,  Br  =  Bh,  we  find  Tr  =  — =  0.685Th  and 
Pm  =0381. 


The  mean  values  and  the  (cross-)variances  of  the  look  images  are  given  by 
lu  =Il  fori  =  1,2, . Ll 


is,  =is=iL/p„  for  j  =  1,2 . Ls 

The  (cross-)variances  are 

(II,  -iL,)(lLj  =  =piL,ii.j  var(lL) 


(II,  -iLi)(ls,-isj)-Pi^i3  ai^o,  _p,^  , 

r  m  Hm 

(is,  -is.)(ls,  -is, )  =  Pi,, I, 


The  multi-look  intensity  image  is 

Ll  Ls 

Im  =0CL]^lLi  +CCs^Isj 

i=l  j=l 
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Thus,  the  mean  value 


-  _  *-’5  _  _ 

Im  =aL^lL.+as^  Isj=aL^iLi+as^isj=aLLLiL+asLs 

i=l 


j=l 


J=« 


P„ 


This  mean  value  is  chosen  equal  to  the  mean  value  of  a  large  look. 
Thus 

oClLl^I-,  “F ocgLg  "  ’  ”  II 
Pm 

Hence 

a,LL+^  =  l 
Pm 

The  variance  is  obtained  by 


var(lM)  =  (lM-lM)  =|ccLX(^Li  -lLi)  +  as^(lsj-Isj)| 

Ll  Ll  Ll  Ls  Ls  Ls 

=  +2aLasX£ar,,o,,.p,,^,,.  +as^£X®'si<^>siPis.isj 

i=l  j=l  i=l  j=I  i=l  j=l 

=  var(l,)  +^«LasXXp>.>.,  +f-ZZp.s,^ 

i=l  j=i  Hm  i=l  j=i  Pm  i=i  j=i 


Hence 
ENL  =  7 


2  Lj  Lj; 


p  “L^sX^PlL^Isj  p2  SXpIsilsj 

i=l  j=i  Pm  J  Pm  i=i  j=i  ^ 


3.4.2  Ll  overlapping  large  Hamming  looks  and  2  small  rectangular  looks 

Figure  16  shows  the  positions  of  the  various  looks  associated  with  this  technique, 
referred  to  as  the  HL7RS2  technique.  Comparing  with  figure  12  in  section  3.3.2 
we  see  that  the  small  Hamming  weighted  looks  are  replaced  by  small  rectangular 
ones.  The  expectation  is  that  by  using  rectangular  weighting  functions  for  the 
extreme  looks  rather  than  Hamming  weighted  ones  we  are  able  to  improve  the 
radiometric  resolution  via  the  edge  looks  without  having  to  sacrifice  geometrical 
resolution. 

In  section  3.3.2  we  used  edge  looks  with  half  the  look  size  of  the  large  ones.  The 
resolution  of  a  rectangular  look  with  this  look  size  is  better  than  the  Hamming  look 
resolution  but  still  worse  than  the  resolution  associated  with  the  large  Hamming 
look.  To  obtain  the  same  geometrical  resolution  we  have  to  use  a  larger  rectangular 
look  size.  Below  we  consider  two  cases:  one  case  where  the  resolution  of  an  edge 
look  equals  that  of  a  large  Hamming  look  and  one  where  the  look  size  of  the  edge 
look  is  one  half  of  that  of  a  large  look. 
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Note  that  by  choosing  the  3  dB  resolution  of  the  rectangular  edge  look  equal  to  the 
resolution  of  the  large  Hamming  look  the  sidelobe  performance  of  the  former  is 
worse. 


large  Hamming  looks 


)' - a-*- - 1 - ^ - i - 1 _ i _ \ _ I 

0  0.5  1  1.5  2  2.5  3  3.5  4 
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Figure  16:  Look  configuration  of  the  HL7RS2  multi-look  technique. 
As  derived  in  section  3.4. 1  we  have 

Im  =  II 


a 


lI-l 


ctsLc  , 
-l--r  =1 


ENL  = 


1 


L,  L, 


i=l  j=I 


+^££p|,<s, 

i=l  j=l 


r  2  ^ 


m  j=i  j=i 


and 


p„= 0.39745^ 


As  an  example,  we  take  an  overlap  of  50%  for  the  large  Hamming  windows,  the 
number  of  independent  looks  Lj  =  4,  Ll  =  2Li-l  and  Ls  =  2. 

The  figures  17  through  18a,b  show  the  results  where  the  resolution  of  the 
rectangular  edge  look  is  equal  to  the  resolution  of  the  large  Hamming  look.  As 
derived  in  section  3.4.1  then  Tr  =  0.685  Th  and  Pm  =  0,581. 


Figure  17  shows  the  ENL  and  3  dB  values  as  a  function  of  as- 
The  maximum  ENL  value  occurs  for  as  =  0.0351  and  is  equal  to  6.70.  This  value 
is  less  than  7.6,  the  value  obtained  with  Hamming  edge  looks.  However,  in  the 
latter  case  the  resolution  is  less,  1.38  compared  to  1.31. 

Compared  to  the  IML  technique  we  notice  that  the  HL7RS2  technique  offers  a 
slight  advantage  in  the  low  as  range.  The  ENL  values  are  almost  equal  but  no 
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resolution  is  sacrificed.  Also,  the  efficiency  is  higher,  only  9  looks  have  to  be 
computed  compared  with  22  in  the  IML  case. 


Ll=7,  50%  overlap  large  Hamming  looks,  Ls=2  small  rectangular  edge  looks,  bm=0.581 
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Figure  17:  The  Equivalent  Number  of  Looks  and  the  3  dB  resolution  as  a  function  ofag.  All  looks  same 
resolution,  =  0.581. 

Figures  18a  and  b  show  the  azimuth  impulse  response  of  the  HL7RS2  technique 
for  three  values  of  ag.  The  value  as  =  0  yields  the  result  of  a  single  look  Hamming 
weighted  image  with  Th  =  Ta/4.  The  other  extreme  value  (ul  =  0)  yields  the  result 
of  a  single  look  rectangular  weighted  image  with  Tr  =  0.685Ta/4.  The  bottom 
curve  shows  the  result  for  as  =  0.035 1  which  corresponds  with  the  maximum  ENL 
value  attainable.  Note  that  although  the  3  dB  resolution  is  equal  to  the  resolution 
of  the  all  Hamming  case,  the  sidelobe  character  is  worse  due  to  the  higher 
sidelobes  of  the  rectangular  look. 

Figure  19  shows  the  result  obtained  with  the  HL7RS2  technique  for  Tr  =  Th/2, 

Pm  =  0.795.  This  figure  can  be  compared  with  Figure  13,  results  from  the  HL7HS2 
technique.  In  both  figures  we  observe  the  as  dependence  of  the  resolution  due  to 
the  inequality  of  the  resolution  of  small  and  large  looks.  The  HL7RS2  technique 
yields  a  maximum  ENL  value  of  7.52  for  as  =  1.356.  The  values  for  the  HL7HS2 
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technique  are  ENL  =  7.56  for  ag  =  0.2.  Thus  the  HL7RS2  technique  yields  a 
slightly  smaller  ENL  value.  However,  the  resolution  is  better,  3  dB  value  is  1.36 
compared  with  1.38,  the  HL7HS2  value. 


Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 


time  arbitrary  unit  ->  aL  =  0.1429  aS  =  0  3-dB  BW  =  1.31  /TL 


time  arbitrary  unit  ->  aL  =  0  aS  =  0.2927  3-dB  BW  =1.31/  TL 


time  arbitrary  unit  ->  aL  =  0.1257  aS  =  0.0351  3-dB  BW  =  1.31  /  TL 


Figure  18a:  Azimuth  impulse  response  of  the  HL7RS2  multi-look  technique  for  several  values  of  ag,  logarithmic 
scale. 
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Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 


Figure  18b:  Azimuth  impulse  response  of  the  HL7RS2  multi-look  technique  for  several  values  of  ag,  linear  scale. 
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3.5  Visual  illustrations  of  multi-look  techniques 

The  results  of  various  multi-look  techniques  are  shown  in  the  form  of  the  images 
R2,  R3  and  R4,  The  techniques  have  been  implemented  in  the  Generic  SAR 
Processor  (GSP)  of  TNO-FEL.  The  GSP  is  a  software  package  that  runs  on  a  SUN 
workstation  and  is  able  to  process  both  airborne  and  satellite  SAR  data.  The  raw 
SAR  input  data  is  part  of  raw  image  data  of  the  ERSl  satellite  SAR.  The  scene 
consists  of  the  entrance  to  the  port  of  Rotterdam  and  part  of  its  environment,  A 
map  of  the  area  is  shown  by  Rl. 

Output  of  the  processor  are  7  single-look  images  (LO,  LI,  L2,  L3,  L4,  L5  and  L6). 
The  images  have  50%  overlap.  Hamming  weighting  is  used.  The  slant  range 
resolution  is  9  m  and  the  azimuth  resolution  is  32  m. 

Image  R2  shows  the  central  single-look  image  L3.  The  speckle  noise  is  clearly 
visible,  masking  weak  details  and  interrupting  line  structures.  The  interruption 
occurs  as  a  low  speckle  value  is  realised. 

Image  R3  shows  the  multi-look  image  consisting  of  the  combination  of  the  4  non 
overlapping  looks:  LO,  L2,  L4  and  L6.  The  corresponding  equivalent  number  of 
looks  ENL  =  4.  The  speckle  reduction,  compared  to  R2,  is  clearly  visible,  more 
weak  details  emerge  and  very  few  line  structures  are  interrupted.  The  probability 
of  occurrence  of  a  low  speckle  value  in  each  of  the  four  images  is  low. 

Image  R4  shows  the  multi-look  image  consisting  of  the  combination  of  all  7 
single-look  images,  LO  through  L6.  The  corresponding  equivalent  number  of  looks 
is  ENL  =  6.4.  Again  an  improvement  can  be  observed  albeit  less  than  going  from 
R2  to  R3. 

Some  remarks  are  in  order.  First,  visual  quality  is  higher  on  a  monitor,  some  is  lost 
by  the  copying  process.  Second,  speckle  reduction  is  advantageous  to  both  men 
and  machine.  A  human  interpreter  sees  more  weak  details.  A  computer  algorithm 
can  use  the  improved  signal-to-noise  ratio  to  perform  better.  Degrees  of  speckle 
reduction  insignificant  to  man  can  still  be  useful  for  very  demanding  applications. 
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R2:  Single  look  oflA  independent  looks. 
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Multi-look  image,  combination  of  4  independent  looks. 
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Multi-look  image,  combination  of?  overlapping  looks. 
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4.  The  azimuth  impulse  response 


In  order  to  assess  a  multi-look  technique  two  values  are  important:  the  Equivalent 
Number  of  Looks  (ENL)  and  the  resolution  measured  by  the  3  dB  width  of  the 
main  lobe  (B  W)  of  the  azimuth  impulse  response.  The  impulse  response  of  the 
multi-look  technique  is  a  weighted  sum  of  the  responses  of  the  large  and  small 
looks. 

i=l  j-1 

The  main  lobe  width  BW  is  determined  via  an  iterative  procedure.  For  an  accurate 
result  a  densely  sampled  impulse  response  Im  is  required.  If  the  constituent 
response  sequences  II  and  Is  are  given  in  formula  form,  the  samples  are  readily 
computed.  In  this  report  only  the  rectangular  and  Hamming  weighting  functions 
are  considered  and  the  corresponding  response  functions  are  analytically  known,  at 
least  the  main  lobe  part.  So,  we  have  used  the  known  formulae.  If  these  formulae 
were  not  known,  a  numerically  obtained  response  sequence  must  be  used.  In  this 
case  the  sequence  must  be  heavily  interpolated  to  yield  a  result  which  can  be  used 
for  a  sufficiently  accurate  BW  value. 

In  the  report  a  number  of  azimuth  impulse  responses  are  shown.  The  main  lobe  of 
these  curves  were  used  to  determine  the  BW  value.  Use  of  the  analytical 
description  of  the  main  lobe  is  correct,  however,  the  description  is  not  correct  for 
the  sidelobes  [ref.  5,  section  7.6].  The  derived  formulae  are  based  on  the 
assumption  that  the  linear  FM  spectrum  has  a  rectangular  amplitude  distribution. 
The  amplitude  ripple  of  the  actual  spectrum  can  be  related  to  the  existence  of 
distortion  sidelobes  that  can  be  considerably  greater  than  the  -40  dB  level  that  can 
be  obtained  by  the  Hamming  weighting  (see  Figure  7.22  in  ref.  5). 

The  Figures  20a  and  20b  show  the  azimuth  impulse  responses  of  the  HL7HS2 
multi-look  technique.  These  curves  were  obtained  by  computer  simulation 
(digitally  computed  cross  correlation  sequence).  Comparing  these  figures  with  the 
ones  shown  by  the  figures  15a  and  b,  which  are  based  on  the  analytical  description 
of  the  main  lobe,  we  see  a  marked  difference  in  sidelobe  behaviour. 

The  computer  simulated  curves  were  obtained  by  using  the  ERS-1  satellite 
parameters.  The  full  resolution  azimuth  response  uses  a  bandwidth  of 
Ba  =  1300  Hz  and  an  integration  time  of  Ta  =  0.6  sec,  yielding  a  bandwidth  time 
product  of  about  800. 

The  spectrum  corresponding  to  such  a  high  BT  product  shows  fairly  well  a 
rectangular  amplitude  distribution  and  a  very  good  sidelobe  behaviour  (not  shown 
here).  However,  the  large  and  small  looks  have  BT  products  which  are  much 
lower.  For  instance,  the  small  look  spans  Ta/8  seconds.  The  azimuth  frequency 

rate  is  —  thus,  the  look  has  a  BT  product  of  only  =  =  12  .5 

Tjj  8  Ta  8  64 

Therefore,  the  spectrum  is  far  from  rectangular  and  a  decreased  sidelobe  character 
can  be  expected. 
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Thus,  we  conclude  that  the  shown  azimuth  impulse  responses  must  be  interpreted 
with  care,  the  main  lobe  part  is  correct,  however,  the  sidelobe  behaviour  is  too 
optimistic. 


Azimuth  impulse  response  of  multilook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 


amplitude  amplitude  amplitude 
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Azimuth  impulse  response  of  multiiook  l=aL*sum  of  7  looks  IL  +  aS*sum  of  2  looks  IS 
0.4 
0.2 
0 

50  100  150  200  250  300  350  400  450  500 

time  arbitrary  unit  ->  aL  =  0.1429  aS  =  0  3-dB  BW  =  1 .3073  /  TL 


time  arbitrary  unit  ->  aL  =  0.1 143  aS  =  0.2  3-dB  BW  =  1 .3793  /  TL 


Figure  20b:  Azimuth  impulse  response  of  the  HL7RS2  multblook  technique  for  several  values  of  ag,  linear  scale. 
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5.  Conclusions 

We  have  evaluated  the  effects  on  SAR  imagery  of  a  number  of  different  speckle 
reduction  techniques.  Each  technique  uses  a  number  of  single  look  intensity 
images  of  which  the  weighted  sum,  the  multi-look  image  is  formed.  The 
constituent  single  looks  differ  in  size  T  (B)  and  weighting  function  used.  It  is 
shown  that  the  speckle  behaviour  of  each  look  is  the  same  and  that  they  differ  in 
geometric  resolution,  i.e.  the  width  of  the  main  lobe  of  the  impulse  response. 

Since,  the  lower  geometric  resolution  looks  do  not  make  efficient  use  of  the  signal 
energy,  speckle  reduction  can  be  obtained,  increase  in  the  ENL’s,  by  multi-look 
technique.  By  using  overlapping  looks,  a  more  efficient  use  is  made  of  the 
available  bandwidth. 

We  have  evaluated  three  types  of  techniques:  techniques  using  equal  looks  with  a 
rectangular  weighting  function,  techniques  that  use  equal  looks  with  a  Hamming 
weighting  function,  and  techniques  that  use  two  kinds  of  weighting  function  sizes. 
It  has  been  shown  that  the  last  techniques  are  able  to  fill  the  left  and  right 
extremities  of  the  signal  bandwidth  by  positioning  there  small  looks.  However,  the 
general  pattern,  trading  geometric  resolution  for  radiometric  resolution,  holds, 
since  the  smaller  looks  have  a  lower  geometric  resolution.  The  IML  technique 
forms  the  multi-look  image  by  computing  the  weighted  sum  of  the  large  and  small 
looks.  By  vaiying  the  factor  as  continuously  a  continuous  trading  of  geometric  and 
radiometric  resolution  is  obtained. 

As  a  measure  of  geometric  resolution,  the  3-dB  width  of  the  main  lobe  of  the 
impulse  response  is  used.  However,  two  responses  can  have  the  same  3-dB  value 
and  differ  in  detailed  main  lobe  and  sidelobe  structures.  The  effects  of  different 
impulse  response  shapes  on  image  interpretation  may  be  important. 

We  have  analysed  the  techniques  assuming  a  flat  azimuth  spectrum.  The  effects  of 
the  shape  of  the  spectrum  can  be  analysed  for  each  individual  system  and  a  more 
realistic  value  of  the  ENL  can  be  calculated.  A  potential  drawback  of  the  multi¬ 
look  techniques  yielding  improved  ENL  values  is  that  many  more  looks  have  to  be 
computed.  However,  there  may  be  cases  where  the  utmost  attainable  ENL  value  is 
a  necessity. 
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Appendix  A  The  correlation  coefficients  pxy  and  pj^j^ 

A.l  The  correlation  coefficient  pxy  of  two  complex  valued  single  look 
images 


Pc 


exp(j  phi(ta)) 


-Ta/i 


Available  integration  time  Ta 


II7 - 

Ta/2 


w(ta) 


->  ta 


The  cross  spectral  density  of  two  complex  single  look  images  is 
P.,(f)=H,(f)H;(f)Pe(f) 

where  Hx(0  and  Hy(f)  are  the  frequency  transfer  functions  of  the  x  and  y  outputs 
respectively,  and  Pc(f)  is  the  azimuth  clutter  power  density  spectrum. 

The  cross-covariance  Pxy  -  ~  x)(y  -  y)  is  given  by 

Pxy  =  Rxy(0)  =  j  Pxy(f)df  =  Pj  H,  (f)H;(f)df  =  P,  J  hx(t)h;(t)dt 
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The  azimuth  reference  function  is 

h(t)  =  ej‘P(')  -Ta/2<t<T,/2 

Hence  the  look  reference  functions  are 

hx(t)  =  Wa(t-tcJh(t)  t,,  -Tx/2<t<te,  +Tx/2 

(0  ~  ''^y(^  ~  (0  ^Cy  ~  Ty/^  ^  t  ^  +  Ty/2 

where  Wx( )  and  Wy( )  are  the  (real)  window  functions. 

Substitution  gives 

^xy  =PcJwx(t-tcx)h(t)wy(t-tcy)h*(t)dt 

M'Xy  “  Pc  j* '^x(t  tcx)Wy^t  —  tcy  ^dt 

Without  loss  of  generality  we  take  Tx  >  Ty  and  tx  =  0. 

Thus 

l^xy  ~  Pc  J*  ^x  (O'^y  ~  ^Cy  ) dt 

The  variance  of  x  and  y  are 

var(x)  =  (x  -  x)(x  -  x)  =PcJwx^(t)dt 

var(y)  =  (y-y)(y-y)  =PcJw)?(t)dt 
Hence,  the  correlation  coefficient  p^y  =  is 

OxOy 

jwx(t)Wy(t-tcy)dt 

■Jjwx^(t)dt,JJw/(t)dt 

A.2  The  correlation  coefficient  pi^i^  of  two  intensity  images 

The  statistics  of  the  real  and  imaginary  part  of  the  complex  valued  images  x  and  y 
are 

Xr  =  Xj  =0 


XfXi  =0 
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var(xr)  =  var(xi)  =  ax^ 
y^  =  Ti=0 

y^=o 

var(yr)  =  var(yi)  =  a/ 


Xryr  =Xiyi 


Xryi  =-Xiyr 

The  stochastic  variables  Xj,  yj-  and  y\  have  a  Gaussian  probability  density 
function.  Thus  the  intensity  has  an  exponential  pdf  with 

h=2a^ 

var(Ix)  =  Ix'=40x‘' 

Thus,  the  statistics  of  the  complex  valued  images  are 
X  =  Xr  +  jXi  =0 

var(x)  =  |x  -  x^^x  ”  xj  =  xx*  =  x^^  +  x^^  =  2o^ 
y  =  yr+jyi  =0 

var(y)  =  (y  -  y)(y  -  y)  =  yy*  =  +  y-^  =  2<Jy 

l^xy  =(x-x)(y-y)  =xy*  =  Xryr+Xiyi+j(xiyr-x^)  =  2(xryr+jXiyr) 


The  correlation  coefficient  between  x  and  y  is  given  by 
^var(x)^var(y)  0,0, 


Hence 


- 2 

Xryr 


+  Xiyr^ 


^  2^  2 
Ox  Oy 


The  real  valued  intensity  images  are 

Ix  =Xr^+Xi^ 


A.3 
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iy=yr'+yi^ 

Thus,  the  statistics  of  the  intensity  images  are 

Ix  =  Xr^  +  Xj^  =  2a^ 

ly='^+77  =  2ay^ 

var(Ix)  =  (lx  -Ix)^  =  4ax‘‘ 

var(l,)  =  (ly-Iyf  =  4a; 

The  cross-covariance  is 

“  (ix  “  “ly)  =  Ixly  ““ixly  =  Ixly  ~4ax^C7y^ 

Ijy  =  (xr^  +  Xi^)(yr^  +  yi^)  =  Xr^yr^  +  Xr^i^  +  +  XjVi^ 

Since  and  yr  are  Gaussian 

Xr  Yr  =Xr  +2Xryr 

Substitution  yields 
IxIy“4^Xryr  +Xiyr  j 
Thus 


Hence 

PMy  =|pxy| 

The  correlation  coefficient  of  the  intensity  images  is  equal  to  the  modules  squared 
of  the  correlation  coefficient  of  the  complex  images. 
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